Attractor Solution in Coupled Yang-Mills Field Dark Energy Models 
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We investigate the attractor solution in the coupled Yang-Mills held dark energy models with 
the general interaction term, and obtain the constraint equations for the interaction if the attractor 
solution exists. The research also shows that, if the attractor solution exists, the equation-of-state of 
the dark energy must evolve from w y > to w y < — 1, which is slightly suggested by the observation. 
At the same time, the total equation-of-state in the attractor solution is Wtot = — 1, the universe is a 
de Sitter expansion, and the cosmic big rip is naturally avoided. These features are all independent 
of the interacting forms. 



I. INTRODUCTION 

Dark energy problem has been one of the most ac- 
tive fields in the model cosmology^ince the discovery of 
accelerated expansion of universe [H 0, S|- I n the obser- 
vational cosmology, the equation-of-state (EOS) of the 
dark energy Wde = Pde/Pde plays a central role, where 
Pde and are its pressure and energy density, respec- 
tively. To accelerate the expansion, the EOS of dark 
energy must satisfy Wde < —1/3. The simplest candi- 
date of dark energy is a tiny positive time-independent 
cosmological constant A, whose EOS is — 1. However, 
it is difficult to understand why the cosmological con- 
stant is about 120 orders of magnitude smaller that its 
natural expectation, i.e. the Planck energy scale den- 
sity. This is the so-called fine-tuning problem. Another 
puzzle of the dark energy is the first cosmological coin- 
cidence problem Q, namely, why does our universe be- 
gin the accelerated expansion recently? why are we living 
in an epoch in which the dark energy density and the 
dust matter energy density are comparable? This prob- 
lem becomes very serious especially for the cosmological 
constant as the dark energy candidate. The cosmologi- 
cal constant remains unchanged while the energy densi- 
ties of dust matter and radiation decrease rapidly with 
the expansion of our universe. Thus, it is necessary to 
make some fine-tuning. In order to give a reasonable 
interpretation to the first cosmological coincidence prob- 
lem, many dynamical dark energy models have been pro- 
posed as alternatives to the cosmological constant, such 
as quintessence 0, phantom Q, k-essenceQ, quintomQ 
etc. 

Recently, by fitting the SNe la data, marginal evidence 
for Wde(z) < —1 at redshift z < 0.2 has been found. In 
addition, many best fits of the present values of Wd e are 
less than —1 in various data fittings with different pa- 
rameterizations. The present observational data seem to 
slightly favor an evolving dark energy with Wd e crossing 
— 1 from above to below in the near pastQ. In has been 
found that the EOS of dark energy Wde cannot cross the 
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so-called phantom divided Wde = — 1 for quintessence, 
phantom or k-essence alone [loj. A number of works have 
discussed the quintom models [8|, which is an combina- 
tion of a quintessence and a phantom. Although many 
of these models provide the possibility that Wd e can cross 
— 1, they do not answer another question, namely, why 
crossing phantom divided occurs recently? Since in many 
existing models whose EOS can cross the phantom di- 
vide, Wde undulated around —1 randomly, why are we 
living in an epoch Wd e < — 1? It is regarded as the sec- 
ond cosmological coincidence problem 11 1 . 

As well known, the most frequently used approach to 
alleviate the first cosmological coincidence problem is the 
tracker field dark energy scenario [12j. The dark energy 
can track the evolution of the background matter in the 
early stage, and only recently, the dark energy has nega- 
tive pressure, and becomes dominant . Thus the current 
condition of the dark energy is nearly independent of the 
initial condition. If the possible interaction between the 
dark energy and background matter [l3j is considered, the 
whole system (including the background matter and dark 
energy) may be eventually attracted into the scaling at- 
tractor, a balance achieved, thanks to the interaction. In 
the scaling attractor, the effective densities of dark en- 
ergy and background matter decrease in the same man- 
ner with the expansion of our universe, and the ratio 
of dark energy and background matter becomes a con- 
stant. So, it is not strange that we are living in an epoch 
when the densities of dark energy and matter are com- 
parable. In this sense, the first cosmological coincidence 
problem is alleviated. On the other hand, if the scal- 
ing attractor also has the property that its EOS of dark 
energy is smaller than —1, the second cosmological coin- 
cidence problem, if existing, is also alleviated at the same 
timefll]. However, this is impossible in the interacting 
quintessence or phantom scenario. 

Recently, a number of authors have discussed another 
class of models, which are based on the conjecture that a 
vector field can be the origin of the dark energy [hH H5|. 
and have different features to those of scalar field. In the 
Refs.0,[TiJ[H[IIj], it is suggested that the Yang-Mills 
(YM) field can be a kind of candidate for such a vector 
field. Compared with the scalar field, the YM field is 
the indispensable cornerstone to particle physics and the 
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gauge bosons have been observed. There is no room for 
adjusting the form of effective YM lagrangian as it is pre- 
dicted by quantum corrections according to field theory. 
In the previous works [H, [l^, |2(j, we have investigated 
the 1-loop YM field case and found attractive features: 
The YM field dark energy models can naturally realize 
the EOS of w y > —1 and w y < —1, and the current state 
of the YM dark energy is independent of the choice of 
the initial condition. The cosmic big rip is also avoided 
in the models. 

In the recent works pll . |22| , the 2-loop and 3-loop YM 
field dark energy are also considered. Although these 
cases are much more complicated than the 1-loop case, 
they have not brought new feature for the evolution of 
the universe. So in this work, we shall only focus on the 
YM field with 1-loop case. 

In this work, the cosmological evolution of the YM 
dark energy interacting with background prefect fluid is 
investigated. In fact, gauge fields play a very important 
role in, and are the indispensable cornerstone to, parti- 
cle physics. All known fundamental interactions between 
particles are mediated through gauge bosons. Generally 
speaking, as a gauge field, the YM field under considera- 
tion may have interactions with other species of particles 
in the universe. However, unlike those well known inter- 
action in QED, QCD, and the electron- weak unification, 
here at the moment we do not yet have a model for the 
details of microscopic interactions between the YM field 
and other particles. In this work, instead of considering 
some specific assumed interactions between YM field and 
matter and radiation, which has adopted in [2^, [U, H2| , 
we shall consider YM dark energy model with a general 
interacting term, and investigate the general feature of 
the attractor solution. 

This paper is organized as follows. In Sec. 2, we give 
out the equations of the dynamical system of the inter- 
acting YM field dark energy models, and discuss the gen- 
eral features of the interacting models. In Sec. 3, we con- 
sider three special cases of the interaction terms and the 
holographic YM dark energy models, and investigate the 
constraints of these interaction terms. Finally, brief con- 
clusion and discussion are given in Sec. 4. 



II. DYNAMIC SYSTEM OF INTERACTING 
YANG-MILLS DARK ENERGY 



attractive features of this effective YM lagrangian include 
the gauge invariance, the Lorentz invariance, the correct 
trace anomaly, and the asymptotic freedom[24j]. With 
the logarithmic dependence on the field strength, ^ e ff 
has a form similar to he Coleman- Weinberg scalar effec- 
tive potential 12711. and the Parker-Raval effective grav- 
ity lagrangian[28]. The effective YM field was firstly 
put into the expanding Robertson- Walker (R-W) space- 
time to study inflationary expansion[16] and the dark 
energy [l7|. We work in a spatially flat R-W spacetime 
with a metric 



ds 2 = a 2 (r)(dr 2 - 8ijdx l dx j ), 



(2) 



where r = J(ao/a)dt is the conformal time. For simplic- 
ity we study the SU(2) group and consider the electric 
case with B 2 = 0. The energy density and pressure of 
the YM field are given by 



Py 



F 2 

T (e + 6). 



E 2 

p » = T U 



where the dielectric constant is given by 
e = bin - 



and the EOS is 



U'u 



Py_ = y-t 

Py 3y + 3 ' 



(3) 



(4) 



(5) 



where y = e/b = ln|^-|. At the critical point with the 
order parameter F — k 2 , one has y = and w y = — 1, 
the universe is in exact de Sitter expansion (l6| . Around 
this critical point, F < n 2 gives y < and w y < — 1, 
and F > k 2 gives y > and w y > — 1. So in the YM 
field model, EOS of w y > —1 and w y < — 1 all can be 
naturally realized. When y ^S> 1, the YM field has a 



state of i 



1/3, becoming a radiation component. The 



effective YM equations are 

<9 M (a 4 e F a ^) + / abc /l^(a 4 e F c ^) = 0, 



(6) 



the v = component of which is an identity, and the 
v = 1,2,3 spatial components of which reduce to 



d T (a 2 eE) = 0. 



(7) 



The effective YM field cosmic model has been discussed 
m Refs.[H,[il [l8l[l9l . The effective lagrangian up to 1- 
loop order is 24 T25 1 



£eff = 7; Fin 



F 



en' 



(1) 



where b — \1N/2Att 2 for the generic gauge group 
SU(N) is the Callan-Symanzik coefficient [2y|. F — 
— {l/2)F® v F a,lu plays the role of the order parameter of 
the YM field, k is the renormalization scale with the di- 
mension of squared mass, the only model parameter. The 



In this work we will generalize the original YM dark 
energy model to include the interaction between the YM 
dark energy and dust matter. We assume the YM dark 
energy and background matter interact through an inter- 
action term Q, according to 



Py + 3H{p y +p y ) = -Q, 



Pv 



3Hp m — Q, 



(8) 



(9) 



which preserves the total energy conservation equation 
Ptot + 3H(ptat + Ptot) — 0. It is worth noting that the 
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equation of motion should be changed when Q ^ 0. 
We introduce the following dimcnsionless variables 



2p n 



6k- 



2Q 



(10) 



where / is the function of x and y. By the help of the 
definition of y, the evolution equations © and © can 
be rewritten as a dynamical system, i.e. 



4y f(x,y) 



2 + y (2 + y)eV 
x' = -3x + f(x,y), 



(11) 
(12) 



here, a prime denotes derivative with respect to the so- 
called e-folding time N = In a. The fractional energy 
densities of dark energy and background matter are given 
by 



(i + y)e y 
(1 + y)ev + x' 



and Qrn = 



(1 + y)e* +x' 



(13) 



We can obtain the critical point (y c , x c ) of the au- 
tonomous system by imposing the conditions y' c — x' c — 
0. From the equations pip and (|12p . we obtain that the 
critical state satisfies the following simple relations 



3x c — f(x c , ?/c): 



3x c 



-4y c e y 



(14) 



(15) 



so we can get the critical state (y c ,x c ) by solving these 
two equations. In order to study the stability of the criti- 
cal point, we substitude linear perturbations y — > y c + 5y 
and x — > x c + Sx about the critical point into dynamical 
system equations (jlll) and (TT2")) and linearize them, and 
obtain two independent evolutive equations, i.e. 



(Sy 1 
\Sx' 

where 



M 



Si) 
5x 



Gy -I - Ry R X 

fy fx ~ 3 



5x 



R y = dR/dy\ {y=yr X=Xa) , 



(16) 



and the definitions of R x , f y , f K and G y are similar. The 
functions G and R are defined by 



G = G(y) = 



4y 
2 + y' 



R = R(x,y) 



f(x,y) 
(2 + y)ev ' 



which are used for the simplification of the notation. The 
two eigenvalues of the coefficient matrix M determine the 
stability of the corresponding critical point. The critical 
point is an attractor solution, which is stable, only if both 
the these two eigenvalues are negative (stable node), or 



real parts of these two eigenvalues are negative and the 
determinant of the matrix M is negative (stable spiral) , 
which requires that the critical point satisfies the follow- 
ing inequalities 



Gy + Ry + f x - 3 < 0, 



(17) 



[R*fy ~ (/x - 3)(G y + Ry)] [(G y +i? y -/ x +3) 2 -4i? x / y ] < 0. 

(18) 

or it satisfies 



Gy + Ry + f X - 3 < 0, 



0. 



(19) 
(20) 



{G y + R y -f x + 3) 2 - 4i? x / y 

These generate a constraint of the interaction term Q, 
which will be shown in the following section. 

Here we discuss some general features of the attrac- 
tor solutions, regardless the special form of the interac- 
tion term Q. From the expression (fT5|) . we find that 



Substitute this into the formula (TIB")) , one 



obtains 



n (jfe + l)e* = 3 + 3y c 

v (y c + l)e»« + aj c 3-y c ' 

Since < Sl y < 1 , this formula follows a constraint of the 
critical point 

- 1 < y c < 0. (22) 

From the formulae (O and (f2Tj) . we obtain 

QyWy = — 1. (23) 

This relation is kept for all attractor solutions, indepen- 
dent of the special form of the interaction. Since the value 
of fly is not larger than one in the attractor solution, we 
obtain that 



w y < -1, 



(24) 



the EOS of the YM dark energy must be not larger than 
— 1, phantom-like or A- like. Since in the early universe, 
the value of the order parameter of the YM field F is 
much larger than that of k 2 , i. e. y 3> 1, the YM field is 
a kind of radiation component [19l|. However, in the late 
attractor solution, the dark energy is phantom-like or A- 
like. So the phantom divide must be crossed in the former 
case, which is different from the interacting quintessence, 
phantom or k-essence models. 

In order to the investigate the finial fate of the uni- 
verse, we should investigate the total EOS in the uni- 
verse, which is defined by 



Plot Py + Pm n 

Wtot = = j = UyWy, 

Ptot Py + Pm 



(25) 



where p m — is used. From the relation (|23|) . we obtain 
that, in the attractor solution, 



Wfot = -1- 



(26) 



This result is also independent of the special form of the 
interaction. So the universe is an exact de Sitter expan- 
sion, and the cosmic big rip is naturally avoided, although 
the YM field dark energy is phantom-like. 
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III. SEVERAL INTERACTION MODELS 

In the previous section, we find the critical point of the 
dynamical system of interacting YM dark energy mod- 
els satisfies not only the equations in (fT4")l and (fTS")) . but 
also the constraint of (|22|) . It is obvious that the expres- 
sion of (|I4j) depends on the special form of interacting 
term. If the critical point is an attractor, it also satis- 
fies the constraint in (TT7]) and (IB]), or in © and ([2"0")) . 
These relations can give some constraints of the inter- 
action term. In this section, we consider several cases 
with different interaction forms between the YM dark 
energy and background matter, which are taken as the 
most familiar interaction terms extensively considered in 
the literature [II]. 

Case a: Q oc Hp yi which is equivalent to the form 
f(x, y) = a(y + l)e v , where a is a dimensionless constant. 
From the equations (| 14[) and (fT5"|) , we obtain the critical 
point 



Vc = 



(27) 



The constraint in (|2"2")) requires that 



a > 0, (28) 

and the attractor conditions in (|I7p - ([20|) require that 

a > -8. (29) 

So we obtain the constraint of the interaction from, if the 
attractor solution exists, the parameter a satisfies 



a > 0, 



(30) 



and the EOS and the fractional energy density of the YM 
field in the attractor solution are 



w 



I 3 

„ - --( a + 3), Sly - — 



a + 3 



(31) 



respectively. It is obvious that w y < — 1. 

Case b: Q cx H(p y + p m ), which is equivalent to the 
form f(x, y) — f3[(y+l)e v +x], where (3 is a dimensionless 
constant. From the equations (fT4")l and (|I5[) . we obtain 
the critical point 



Vc 



3/3 



4y c 



(32) 



/3-12' 3 
The constraint in (f22|) requires that 

< (3 < 3, (33) 
and the attractor conditions in (|I7 p - (|2"D|) require that 

120 



31 



So the parameter a satisfies 

< (3 < 3, 



(34) 



(35) 



if the critical state is an attractor solution. The EOS 
and the fractional energy density of the YM field in the 
attractor are 



3 _ 3-/? 

Wy ~ /3-3' v ~ 3 ' 



(36) 



respectively, which follows that w y < —1, the YM field 
dark energy is phantom-like or A-likc. 

Case c: Q oc Hp m , which is equivalent to the form 
f(x,y) = jx, where 7 is a dimensionless constant. From 
the equations (fT4]l and (fT5")) , we easily nd that they have 
no solution except that the value of 7 is exactly zero, i.e. 
the case with no interaction. 

Case d: Recently, a number of authors have discussed 
the holographic dark energy, where the holographic prin- 
ciple has been put forward to explain the dark energy. 
According to the holographic principle, the number of 
degrees of freedom of a physical system scales with the 
area of its boundary. In the context, Cohen et al[29l| 
suggested that in quantum field theory a short distant 
cutoff is related to a long distant cutoff due to the limit 
set by formation of a black hole, which results in an upper 
bound on zero-point energy density. In line with this sug- 
gest, Hsu and Li[3(l[3l| argued that this energy density 
could be views as the holographic dark energy satisfying 



M 2 M 2 P L- 



(37) 



where d > is a numerical constant, and Mp = I/v8vrG 
is the reduced Planck mass. L is the size of the current 
universe. Li[3l| proposed that the IR cut-off L should be 
taken as the size of the future event horizon 



. dt 

L = R eh (a) = a / — j- = a 
It a(t) 



da 
lid? 



(38) 



In this letter, we consider the holographic YM field 
dark energy. From the relation (f37|) . one obtains that 



Sin 



Py = f>de = 6MpSl y H 3 



which follows that the interaction form is 



(39) 



f Sh, 



= -2p y Hl^_l-l\-3Hp y (l + w y ), (40) 



where the expression ([8]) is used. This formula is equiv- 
alent to the form 



f{x,y) 



y+ 1 




(y+l) e y. (41) 



From the equations (|I4j) and (fT5"|) , we obtain the critical 
point 



Vc 



3(d 2 - 1) 
3 + d 2 



4y c 



(42) 
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The attractor conditions in (jTTJ) - (J20J) require that 

d < 0, (43) 

which is conflicting with the previous assumption d > 
0. So we get the conclusion: the holographic YM dark 
energy model has no attractor solution. 

IV. CONCLUSION AND DISCUSSION 

In summary, the cosmological evolution of the Yang- 
Mill field dark energy interacting with background mat- 
ter is investigated in this letter. We find the features of 
the interacting YM dark energy models: 

a. The interaction term between the YM dark energy 
model and the matter has a fairly tight constraint, if we 
require that the attractor solution of the model exists. 

b. If the attractor solution exists, the EOS of the YM 
field must evolve from w y > to w y < — 1 or w y — — 1. 

c. The holographic YM dark energy model has no at- 
tractor solution, which is different from other holographic 
models [32f. 

d. In the attractor solution, the total EOS is Wtot = 
— 1, which is independent of the interacting forms. So 
the universe is in a de Sitter expansion, and the cosmic 
big rip does not exist in the models. 

In the interacting YM dark energy models, we should 
notice the "fine-tuning" problem, which is reflected by 
the value of k, the energy scale of the Yang-Mills field 
dark energy models. In the interacting models, the total 
energy density in the universe is 

P to t = ^ = bj ^l(l + y)e y +x}. (44) 



In the attractor solution, we can obtain 

(Hot = h -Y (i - \y)j * v °. (45) 

where the express (15) is used. The value of ptot should 
be not larger than which of the present total energy den- 
sity in the universe [23|, i.e. 

Ptot < 8.099/i 2 x KT^eF 4 , (46) 

which leads to 

/ i \-i/2 
k < 4.18/i x lO- 5 eV 2 (l - -y c J e _yc/2 . (47) 

For a fixed interacting models, where y c can be obtained, 
one can exactly calculate the value of the k, which keeps 
the current energy density of YM dark energy being cur- 
rent observed value. From (47), we find that this energy 
scale k is, as well as the case with free Yang-Mills field 
models, very low compared to the typical energy scales 
in particle physics. 
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